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Abstract 

The lifetime of a system of connected units under some natural assumptions can 
be represented as a random variable Y defined as a weighted lattice polynomial of 
random lifetimes of its components. As such, the concept of a random variable Y 
defined by a weighted lattice polynomial of (lattice-valued) random variables is con- 
sidered in general and in some special cases. The central object of interest is the 
cumulative distribution function of y. In particular, numerous results are obtained 
for lattice polynomials and weighted lattice polynomials in case of independent argu- 
ments and in general. For the general case, the technique consists in considering the 
joint probability generating function of "indicator" variables. A connection is studied 
between Y and order statistics of the set of arguments. 

Keywords: lifetime models, reliability, weighted lattice polynomial, discrete Sugeno inte- 
gral, probability generating function. 



1 Introduction 

In a system consisting of several interconnected units its lifetime can often be expressed in 
terms of those of the components. Assuming that the system and the units are of the crisply 
"on/off" kind, the lifetime of a "serially connected" segment of units (where all units are 
needed for functioning) is the minimum of the individual lifetimes; for a circuit of parallel 
units (one unit is sufficient for work) the lifetime is the maximum of the individual lifetimes. 
In a complex system where both connection types are involved the system's lifetime Y can 
be expressed as a lattice polynomial function (combination of minima and maxima) of the 
component's lifetimes. 
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In addition, there may be "collective upper bounds" on lifetimes of certain subsets of 
units imposed by external conditions (say, physical properties of the assembly) or even 
"collective lower bounds" imposed for instance by back-up blocks with constant lifetimes. 
These upper and lower bounds can be easily modelled by incorporating additional com- 
ponents with constant lifetimes. Indeed, consider a subset with random lifetime X and 
assume that a lower (resp. upper) bound c G M is required on this lifetime, which then 
becomes c V X (resp. c A X). This situation can be modelled by connecting in parallel 
(resp. series) each component of that subset with a component of constant lifetime c. 

Following the above description, the lifetime of such a system involves the lower and 
upper bounds as follows. Denoting by [n] = {1, . . . ,n} the set of indices of the units and 
by Xi the random lifetime of the ith unit {i G [n]), we conclude that 

y = p(Xi,...,x„) 

where p is an n-ary weighted lattice polynomialjl] whose definition is recalled next (see 
Marichal [1]). Let L denote an arbitrary bounded distributive lattice with lattice operations 
A and V and denote by a and h its bottom and top elements. For any a, f3 & L and any 
subset C [n], let e^'^ denote the characteristic vector of S in {a,/5}", that is, the n-tuple 
whose ith component is /3, if 2 G S, and a, otherwise. The special case e^'^ will be denoted 
els- 
Definition 1. The class of weighted lattice polynomials from to L is defined as follows: 

(1) For any k G [n] and any c E L, the projection (xi, . . . , x„) i— > Xk and the constant 
function (xi, . . . , x„,) t— > c are weighted lattice polynomials from L" to L. 

(2) If p and q are weighted lattice polynomials from L" to L, then p A q and p\/ q are 
weighted lattice polynomials from L" to L. 

(3) Every weighted lattice polynomial from L" to L is constructed by finitely many ap- 
plications of the rules (1) and (2). 

It was proved ^ that any weighted lattice polynomial p : L" ^ L can be expressed in 
disjunctive normal form : there exists a set function w : 2^"' L such that 



pix)= V wiS)Af\ 



Xi 



(1) 



Moreover, from among all the possible set functions : 2'"] ^ L fulfilling ([T]), only one is 
nondecreasing. It is defined as 



w{S)=p{e^/), SC[n]. (2) 

The n-ary weighted lattice polynomial defined by a given nondecreasing set function 
-u; ; 2W ^ L will henceforth be denoted p^. 

Remark. Weighted lattice polynomials are strongly related to the concept of discrete 
Sugeno integral [Qj [10]. Indeed, as observed in [1], a function / : L" — L is a discrete 
Sugeno integral if and only if it is a weighted lattice polynomial fulfilling /(e0^) = a and 
/(eXt ) = h. 



^In the rest of the paper, weighted lattice polynomial functions will be called weighted lattice polynomials. 
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Assuming that L C M = [— cx), oo] and that the hfetime Xi of every unit i G [n] is 
L-valued, the hfetime of the system, with possible cohective lower or upper bounds, is then 
given by 

r = p^(Xi,...,x„) 

where pw is a weighted lattice polynomial from to L. 

In those terms, Pr[y > y] is the probability that by the time y the system is still on, 
that is, the system reliability at time y. This probability is in turn determined by the 
pattern of units that are on at the time y, that is, by the indicators of events Xj > y. 

In this paper we investigate the cumulative distribution function (c.d.f.) of Y in case 
of general dependent arguments (see Section 2). Our results generalize both the case of 
independent arguments, which was studied in Marichal [6], and the special situation where 
no constant lifetimes are considered (the weighted lattice polynomial is merely a lattice 
polynomial), which was investigated by Marichal [5] for independent arguments and then 
by Dukhovny [2j for general dependent arguments. 

We also investigate the special case where the weighted lattice polynomial is a symmetric 
function and the case where, at any time, the reliability of any group depends only on the 
number of units in the group (see Section 3). 

2 Cumulative distribution function in the general case 

We assume throughout that L C M, which implies that L is totally ordered. For an L- valued 
random variable X one can introduce a supplementary variable 

Xix) = lnd{X > x) e {0,1}. 

For a set of L- valued random variables Xi, i E [n], we consider a vector of (synchronous) 
indicator variables 

where 

Xi{x) = lnd{Xi > x), i e[n]. 

For any k G {0, 1}", let Sf^ denote the subset of [n] whose characteristic vector in {0, l}'^ 
is K (that is such that es^ = k). 

We have the following fundamental lemma: 

Lemma 2. Let p^ : L"' —>■ L be a weighted lattice polynomial and Y = p^(Xi, . . . ,X„). 
Then 

lnd{Y > y) = lnd{w{S^^y)) > y), y E L. 

Proof. The proof mainly lies on the distributive property of Ind(-) with respect to disjunc- 
tion and conjunction: for any events A and B, 

Ind(AVfi) = Ind(A) Vlnd(S), 
Ind(AA5) = Ind(yl) Alnd(S). 
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By ([T]), we obtain 

lnd{Y>y) = lnd{pUXi,...,X^)>y) 

= V (lnd{w{S)>y)A/\lnd{X,>y) 

SC[n] ^ ies 

= T^ujyixiy)), 

where 71^^^ : {0, 1}" {0,1} is the weighted lattice polynomial defined by the nondecreasing 
set function Uy : 2^ {0, 1} fulfilling 

ujy{S) = lnd{w{S)>y), SC[n]. 

Finally, since x{y) G {0, 1}", by we have 

T^i^yixiy)) = ^yiSx(y)) = lnd{w{S^(^y)) > y). 

□ 

Remark. Interestingly enough, we have observed in the proof of Lemma |2] that the 
identity Y = py^{Xi, . . . , Xn) has its binary counterpart in terms of indicator variables, 
namely 

xviy) = TTu.y{xi{y),---,Xn{y)), y ^ L, 

where Xriy) = Ind(F > y). 

We now state our main result. 

Theorem 3. Let Y = pw{Xi, . . . , X„) and denote its c.d.f. by Fyiy)- Then, for any y & L, 

1-Fy(y)= Yl lnd{w{S,)>y)PT{x{y) = f^), (3) 
Ke{o,i}" 

or equivalently, 

1 - Friy) = E Ind(«;(5) > y) Pr(x(2/) = Ss). (4) 

SC[n] 

Proof. By Lemma [2|, we have 

l-Friy) = E(Ind(y > y)) = E(Ind(«;(%(,)) > y)) 

= ^ E{lnd{w{S,) > y)) Piixiy) = i^) 
Ke{o,i}" 

= Ind(«;(S'^) > y) Piixiy) = f^)- 

kG{0,1}" 

□ 

Also, since 

lnd{wiS^) >y) = l- lnd{w{S,) ^ y), 
and Z]KPr(x(?/) = k) = 1, one can present ([3]) as 

FY{y)= E lnd{w{S,)^y)PT{xiy) = f^). 
Ke{o,i}" 

The formulas of Theorem [3] are quite remarkable. They show that the vector of binary 
status indicators x{y) at the time y is the only thing needed to determine the status of the 
whole system. Moreover, the "dot product" nature of these formulas has, as it will be seen, 
a special significance. 
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Remark. It is noteworthy that an alternative form of formula (jl]), namely 



1 - Friy) = E Ind(«;(,S) > y) Pr V X, ^ y < f\ X, 

SC[n] ^ie[n]\S ies 

was already found through a direct use of the disjunctive normal form ([1]) of pyj (see proof 
of Theorem 1 in Marichal [6] ) . 

In the general case when the random arguments are dependent the probability parts of 
the formulas ([3])-(jl]) can be treated using the technique of probability generating functions 
(p.g.f.). For a random variable X with a c.d.f. F{x) we introduce the p.g.f. of its indicator 
X{x): 

G{z,x) = E(z^(^)) = Ft[x{x) = 0] + zFt[x{x) = 1] 
= F{x) + z[l - F{x)], \z\ ^ 1. 

For a set of (generally dependent) random variables Xi, . . . ,Xn, we denote the joint 
p.g.f. of their indicator variables x(x) by 

G{z,x) = E(f^>)) = e( n^^'^'^'M, ^ 1, « G N- 

^ j=i ^ 

By definition, 

G{z, x)= Yl P^xix) = k) r = Yl Pr(x(x) = es) Y[ z,. (5) 

kG{0,1}" SC[n] ieS 

Since G{z, x) is a multilinear polynomial in zi, . . . , Zn, we immediately obtain 

Pr(x(x) = k) = — ^ G{z, x) 

oz z=6 

which, with ([3]), enables to compute Fy^y) whenever G{z,y) can be provided. 

As a corollary, we immediately retrieve the known expression of Fyiy) in the special 
case when Xi, . . . , Xn are independent (see Marichal [6]). 

Corollary 4. When Xi, . . . ,X„ are independent random variables with c.d.f. 's Fi{x), i G 
[n], we have 

n 

l-Fy{y)= ^^d{w{S^)>y)X{F,{yf-^^[l-F,{y)p. 
Ke{o,i}" «=i 

Proof. By independence, we immediately have 

n 

G{z,y) = Y{[F,{y) + z,{l-F,{y))l 

i=l 

and hence 

""^^ ' l[F,{yY--^[l-F,{y)r. 



z=0 



1=1 



□ 



^Note that, since G(z, x) is a multilinear polynomial, its derivatives can be also obtained by using 
discrete derivatives (see Grabisch et al. [31 §5]). 
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Remark. The assumptions of Corollary H] can be slightly relaxed while keeping a product 
form for Pr{x{x) = k). In fact, assuming that the arguments are dependent but their 
synchronous indicators Xi{x), . . . , Xni^) are independent for any fixed x, we get 

n 

Pr(x(x) = k) = l[FI[x^{x) = 0]'-^^ PT[x^ix) = 1^- 

i=l 

Similarly, in the case when all Xi = Z + Yi, Z and all Yi being independent (random shift), 
the probability density function formula is the one for the independent arguments case 
convoluted with the density function of Z. 

An alternative approach for computing Pr(x(x) = Ss) is given in the following proposi- 
tion, where F{x) denotes the joint c.d.f. of Xi, . . . , X„: 

Proposition 5. For any fixed x L, the set function S i— > Pr(x(x) = Ss) is the Mobius 
transform of the set function S ^ G{es, x) and we have 

G{es,x)=F{ef), SC[n]. 

Proof. Let x E L. By ([5]), we have 

G{es,x) = J2 ^"^(xix) = St), 

TCS 

which establishes the first part of the result. From this latter identity, we obtain 

G{es, x) = ^ X G M \ T and Xi > x \/i e T) 

TCS 

= Pr(Xi ^xWie [n]\S), 

which proves the second part. □ 
The result in Proposition [5] yields an explicit form of G{z, x) in terms of F{eg'^), namely 

G{z,x)= E nef) n^^ n (i-^o- 

SC[n] ieS ie[n]\S 

Moreover, combining Proposition \5\ and Mobius inversion formula, it follows immediately 
that 

Pi{x{x)=es) = Y.i-^f^'^^^G{eT,x) 

TCS 

= E(-l)'"'-'"'^(eT'Y (6) 

TCS 

Remark. We have observed that, since G{z,y) is a multilinear polynomial, it is com- 
pletely determined by the 2" values G{es,y) = F{ef^), S C [n\. Thus, the computation 
of Fyiv) does not require the complete knowledge of the joint c.d.f. F{x), but only of the 
values F(e|'^) for all y G L and all 5 C [n]. 
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3 Selected special cases 

In this final section we consider three special cases that have a natural interpretation in 
terms of eligible applications: 

1. The weighted lattice polynomial is a lattice polynomial (there are no lower or upper 
bounds) . 

2. The weighted lattice polynomial is a symmetric function (the system lifetime is a 
symmetric function of the component lifetimes). 

3. The synchronous indicator variables are exchangeable (at any time, the reliability of 
any group depends only on the number of units in the group). 

For any binary vector k G {0, 1}", we set |k| = Hi- 
3.1 Lattice polynomials 

An important special case arises when the lifetime of a system is obtained as a lattice 
polynomial of units' lifetimes. Recall that lattice polynomials are defined as follows (see 

i): 

Definition 6. The class of lattice polynomials from L" to L is defined as follows: 

(1) For any A; G [n], the projection (xi, . . . , Xn) i— >■ is a lattice polynomial from L" to 
L. 

(2) If p and q are lattice polynomials from L" to L, then p A q and p V q are lattice 
polynomials from L" to L. 

(3) Every lattice polynomial from L" to L is constructed by finitely many applications of 
the rules (1) and (2). 

Thus, a lattice polynomial from L" to L is nothing else than a weighted lattice polyno- 
mial Pu, '■ ^ L with a set function w : 2'"] — {a, b} fulfilling w{0) = a and w([^]) = b 
(see [3]). In this case, for any y E L, we obtain 

Ind(w;(S') > y) = Ind(j3^(e5 '') > y) = Ind(p^(e5 '') = b) Ind(6 > y) 

= TTUes)lnd{b> y), (7) 

where vr^ : {0, 1}" {0, 1} is the (weighted) lattice polynomial defined by the nondecreas- 
ing set function io : 2^"'^ — > {0, 1} fulfilling 

uj{S) = lnd{pui{es^) = b), 5* C [n]. 

Remark. Clearly, vr^ is exactly the weighted lattice polynomial tt^^^ introduced in the 
proof of Lemma [2] computed for the case of lattice polynomials. Moreover, it has the same 
disjunctive normal form as pyj. 

Since Pr(x(6) = k) = except when k = 0, we obtain the following corollary (see 
Dukhovny 0): 
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Corollary 7. Under the assumptions of Theorem\^ and assuming further thatp^ is a lattice 
polynomial, we have 

l-Fy{y)= 7r^(/?) Pr(x(?/) = (8) 

Ke{o,i}" 

in case of generally dependent argument and, in case of independence, 

n 

l-Fy{y)= T^u{K)X{F,{yf-^^\l-F,{yT. 
Ke{o,i}" i=i 

3.2 Symmetric weighted lattice polynomials 

Considering a system with collective upper or lower bounds on subset lifetimes, an impor- 
tant special case arises where the corresponding weighted lattice polynomial is a symmetric 
function, that is, invariant under permutation of its variables. 

Prominent instances of symmetric weighted lattice polynomials are given by the so-called 
order statistic functions, whose definition is recalled next; see for instance Ovchinnikov [7]. 

Definition 8. For any k E [n], the kill order statistic function is the lattice polynomial 
fk : ^ L defined as 

fk{x) = X(fc) = V /\Xi. 

SC[n] ies 
\S\=n-k+l 

(For future needs we also formally define /o = a and fn+i = b.) 

Remark. A system whose underlying weighted lattice polynomial is the order statistic 
function fn-k+i is called a k-out-of-n system] see for instance the book by Rausand and 
H0yland [HI §3.10]. 

The following proposition shows that a weighted lattice polynomial : — > L is 
symmetric if and only if its underlying set function w : 2t"'l {a, 6}, as defined in ([2]), is 
cardinality-based, that is, such that 

|5| = |^'| ^ w{S)=w{S'). 

Equivalently, there exists a nondecreasing function m : {0, 1, . . . ,n} L such that 

w{S) = m{\S\), SC[n]. 

The proposition also provides an expression of an arbitrary symmetric weighted lattice 
polynomial as a "convolution" of L- valued functions on [n]: m{s) and l{s) = fs{x). 

Proposition 9. Consider a weighted lattice polynomial p^ : L" — > L. Then the following 
assertions are equivalent: 

(i) The set function w is cardinality-based. 

(a) There exists a nondecreasing function m : {0, 1, . . . , n} — L such that 

mis) A fn~s+i{x) . 



Pw{x) = V 
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(in) Pui is symmetric. 
Proof, (i) (a) By ([1]), we have 



Pw{x) = y w{S) A /\xi = y m{s) A 

SC\n] ieS ^ s=0 



y /\x, 

L SC[n] ieS . 
\S\=s 



which proves the stated formula, 
(n) =^ (Hi) Trivial. 

{Hi) =^ (i) Follows immediately from ([2]). □ 

Remark. It is noteworthy that Definition [8] and Proposition [9] can also be stated in the 
more general case where L is an arbitrary bounded distributive lattice. 

For symmetric weighted lattice polynomials, Theorem [3] takes the following form. For 

any y E L, set s{y) = min{s, n + 1 : m{s) > y}. 

Theorem 10. Under the assumptions of Theorem and assuming further that the set 
function w is cardinality-has ed, we have 



l-FY{y)= E Pr(|x(2/)l = s). 

s=s{y) 

Proof. When w is cardinality-based, the sum in ([3]) can be rewritten as 

n 

l-Fyiy) = 51 Ind(m(s) > 2/) ^ Pr(x(2/) = 

s=0 kG{0,1}" 
\k\=s 

n 

= Elnd(m(s) > y) Pr(|x(y)| = s). 

s=Q 



(9) 



(10) 
□ 

fll) 



When m{n) ^ y, we have s{y) = n + 1, and ffTU]) reduces to 0. 

Obviously, ([9]) can also be presented as 

s{y)-i 

Fy{y)= Y: Pr{\x{y)\ = s). 

Both Qj and ffTTj) require the knowledge of Pr(|x(y)| = s), which can be calculated 
through the following formula: 

Proposition 11. For any s G {0, 1, . . . , n}, we have 



mxiy)\=s) = l^_^Gitl,y) 



t=o 



Proof. From it follows that 
Gitty) = 



Ke{o,i}" 



E 



?e{o,i}" 



j:prmy)\ = s)f 



s=0 



from which the result follows. 



□ 
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In terms of the joint c.d.f. of Xi, . . . , we also have the following formula: 
Proposition 12. For any s G {0, 1, . . . , n}, we have 



T ) 



:i2) 



Proof. From it follows that 

mm\ = s) = j: Frixiy) = es) = e E(-i)'"""^'^(e"i'') 

5C[n] SC[n] TCS 

\S\=s \S\=s 

= E Fie^') E 

TC[n] 5Dr 
|T|s:s \S\=s 

where the inner sum reduces to (— 1)"*"'^' • □ 

We now derive an expression of PT{\x{y) | = s) is terms of the c.d.f.'s of the order statistic 
functions. 

Consider the random variable X^k) = /fc(Xi, . . . ,X„) and its c.d.f. F(fc)(2/) = Fx^k)iy)- 
It is easy to see that, for the A;th order statistic function fk, we have s{y) = n — k + 1 if 
y < b and s{y) = n + 1 ii y = b. Therefore, from (fTTj) it follows that 

n—k 

Fik)iy) = E Pr(lx(2/)l = s) = Pr(|x(?/)| ^n-k), (13) 

s=0 

from which we immediately derive the following formula: [f] 

Proposition 13. For any s G {0, 1, . . . , n}, we have 

Pr(|x(?/)| = s) = - (14) 

Incidentally, by combining f[T^ and f[T^ . we get an explicit form of F(fc)(2/) in terms of 
F{x); see also David and Nagaraja [1, §5.3]. 



Corollary 14. We have 



FM= E(-i)'''-'=('fLiV( 

5Cfnl \ / 



5C[n 



Proof. Combining f[T^ and ffT5]l . we obtain 



FM = EE E (-irT H"(^"'t") 



s=Ot=OTC[n] ^/ 
|T|=t 

n— fc n— fc 

E E ^(e^) E(-i)^"* 

t=0 TC[n] s=i 
|T|=t 



^In we formally define F(o)(y) = f and i^(„^_i')(j/) = 
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where the inner sum reduces to the following classical binomial identity 



s=0 

Therefore, 



TC[n] 
\T\sin-k 



k-1 

and the result follows from the substitution S = [n]\T. □ 



3.3 Cardinality symmetry of the indicator variables 

A case is now considered where the joint c.d.f. of the arguments Xi, . . . ,X„ is such that, 
for any x G L, Pr(x(x) = k) depends only on that is 

\k\ = \k,'\ =^ Pt{x{x) = k) = Pr(x(z) = k.'). 

Equivalently, for any x E L, there exists a function : {0,1, . . . ,n} ^ [0, 1] such that 

FT{x{x) = K) = g,{\K\), Ke{0,ir. (15) 

In terms of the analogy with system reliability, assumption ffTSj) means that the probability 
that a group of units survives beyond x (that is, the reliability of this group at time 
x) depends only on the number of units in the group, which is why we call assumption 
(1151) "cardinality symmetry". In view of (jH]), this assumption is satisfied, for instance, 
when Xi, . . . ,Xn are exchangeable, that is, when their joint c.d.f. is invariant under any 
permutation of indices. However, f|T5|) itself suggests only that the synchronous indicator 
variables be exchangeable. 

In particular, for independent and identically distributed (i.i.d.) arguments with a 
common c.d.f. F{x), we have (see Corollary H]) : 

Pr(x(x) = k) = F{xY-\% - F{x)f\. 

More generally, one can easily see that when the arguments Xi, . . . , Xn are obtained from 
some i.i.d. random variables li, . . . , y„ by a transformation 

X, = f{Y„Z), ie[n], 

where Z is a vector of some (random) parameters, assumption fll5p still holds true. 

Corollary 15. Under the assumptions of Theorem\^ and assuming further that the joint 
c.d.f. possesses property l[T5\) . we have 

n 

l-Fy{y) = J29yis) E l^diwiS)>y). (16) 

s=0 SC[n] 
\S\=s 

In particular, in case of lattice polynomials, 

n 

1- Fyiy) = Y.9y{s) E ^'^(^)- 
s=0 kG{0,1}" 
\k\=s 
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By Proposition [5], we also have 



i=0 



51 /,^| 



which shows that the set function S (—> Ficg'^) is cardinahty-based. 
For order statistic functions, we have □ 



n-k / \ 

Fik){y) = Y.[A9y{s). (17) 

Indeed, based on (fT^ . we have 

F(fc)(y)=Pr(|x(y)|^n-fc)= 9ym), 

Ke{o,i}" 

|K|^n— A; 

which proves (1171) . 

It turns out that cardinality symmetry is necessary and sufficient for the special relation 
between the c.d.f.'s of Y and order statistic functions provided in Marichal [5], Theorem 8] 
for the case of a lattice polynomial of i.i.d. arguments. Beyond its theoretical value, the 
relation is a simple linear equation and can be used to facilitate computations. 

Theorem 16. Under the assumptions of Theorem\^ the following relations hold true for 
an arbitrary weighted lattice polynomial if and only if the arguments Xi, . . . ,X„, possess 
cardinality symmetry: 

n 

\-Fy{y) = J2^4F(^n-s)iy) - F^n-s+i){y)i (18) 

s=0 

= ^ Ind(w(5) > y), s = 0,...,n. (19) 

\s) SC[n] 
\S\=s 

Proof. (<^=) Assuming (fT5l) . on the strength of (fTTl) . 



n ' 
. s I 



9y{s) = F^n-s){y) - Fi^n-s+i){y), s = 0, ...,n. (20) 



Now, using (|20D in ([TSD, we obtain (118!) -([T9|). 

(^) Assume that f|T8|) and f|T9l) are true for any weighted lattice polynomial and prove 
that f|T5|) must hold true. 

Consider the particular case of a lattice polynomial pw Substituting (JTl) and in 
(HSD-dUD we obtain 

n 1 

^-My) = Y.7n^ E ^^i^) E Pr(x(i/) = /?). (21) 

s=o Ke{o,i}" kg{o,i}" 



''The total probability condition allows to formally extend ([T7]) to the case k = Q. 
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On the other hand, by ([8]) we have 

n 

l-Fy{y) = Y E ^my) = ^) (22) 

s=o Ke{o,i}" 

k|=s 

Therefore, the right-hand sides of 0211) and fl22p coincide for any lattice polynomial Pw In 
particular, they coincide for the lattice polynomial Pw ^ fk- Since 

(tTo; a fk){i^) = 1 -v^ T^ujii^) = 1 and ^ n — k + 1, 

it follows that 

E 7^ E ^^('^) E Pr(x(2/) = '^) 

s=n-fc+l Ke{0,l}" kG{0,1}" 

|k|=s I'?I=* 

n 

= E E TT^l/t) Pr(x(2/) = /i; = 0, ...,n. 

s=n-A;+l kg{0,1}" 
\k\=s 

It implies that for any lattice polynomial it must hold true that 

73- E ^^('^) E = fi) 

\k\=s I'^I=* 

= Y ^U^)'PTixiy) = ^), s = 0,...,n. (23) 

Ke{o,i}" 

|k|=s 

Fix s G {0, . . . ,n} and X C [ri] such that \K\ = s. Applying ( l23il to the lattice polynomial 
7ra;(/t) = /\ieKi^i, it yields that 

"4: E = /«) = Pr(x(2/) = /«') 

VJ '?G{0,1}" 
k|=s 

for any G {0, 1}" such that \k'\ = s, which completes the proof. □ 
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